7 Be(p,7) 8 B astrophysical S-factor from precision cross section measurements 
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We measured the 7 Be(p,7) 8 B cross section from E cm = 186 to 1200 keV, with a statistical-plus- 
systematic precision per point of better than ±5%. All important systematic errors were measured 
including 8 B backscattering losses. We obtain Si7(0) = 22.3 ± 0.7(expt) ± 0.5(theor) eV-b from 
our data at E cm < 300 keV and the theory of Descouvemont and Baye. 

PACS numbers: 26.20+f, 26.65+t, 25.40Lw 



It is now known that electron neutrinos (f e 's) from the 
decay of 8 B in the Sun oscillate into u^'s and/or z^-'s, and 
possibly into sterile v x 's [Q. The v e production rate is 
based on solar-model calculations that incorporate mea- 
sured reaction rates for most of the solar burning steps, 
the most uncertain of which is the 7 Be(p,7) 8 B rate. Im- 
proved production rate predictions are very important 
for limiting the allowed neutrino mixing parameters in- 
cluding possible contributions of sterile neutrinos. The 
astrophysical S-factor Si7(0) for this reaction must be 
known to ±5% in order that its uncertainty not be the 
dominant error in predictions of the solar v e flux Q . 

Si7(0) values based on previous direct measurements 
have quoted uncertainties of typically ±9% or larger ||, 

i iJL @> ! !> ( see also the quoted ±5% results of 
rcf. [|ll| ), while for many of these experiments there are 
unsettled issues such as possible 8 B backscattering losses. 
Indirect Si7(0) determinations based on Coulomb disso- 
ciation and peripheral transfer reactions are also avail- 
able but it is difficult to determine all of their im- 
portant systematic errors. 

We have made a precise determination of Si7(0) using 
a technique that incorporates several improvements over 
traditional methods. We avoided a major difficulty in 
most previous experiments due to uncertain and nonuni- 
form target areal density by using a ^1 mm diameter 
beam magnetically rastered to produce a nearly uniform 
flux over a small ~3 mm diameter target. We directly 
measured the energy loss profile of the target using a nar- 
row 7 Be(a, 7) n C resonance and we determined all impor- 
tant sources of systematic error including the first direct 
measurement of 8 B backscattering losses. 

We used a 106 mCi 7 Be metal target fabricated at TRI- 
UMF and deposited on a molybdenum backing. The 
cross sections were measured using the University of 
Washington FN tandem accelerator with a terminal ion 
source. A proton beam, typically 10 /iA, passed through 
an LN 2 -filled cold trap directly upstream of the target. 
Cryopumps were used for high-vacuum pumping, and 



sorption pumps for roughing. The water-cooled target, 
and a plate with precision-sized circular apertures were 
mounted on opposite ends of a rotating arm. Rotating 
the arm 180° from its horizontal bombardment position 
placed a 3 mm aperture in the beam, and the target ~4.5 
mm from a 450 mm 2 40 micron Si detector that counted 
/3-delayed a's from 8 B decay. In each measurement, the 
arm was rotated through many complete cycles. 

We integrated 3 different beam currents: the current 
striking the target during the bombardment phase, and, 
during the cv-counting phase, the current striking the 
aperture and the current collected in a Faraday cup after 
passing through the aperture. The target arm was biased 
to +300V. The neutral H content of the beam was found 
to be < 10~ 4 , and the cup current changed by < 0.5% for 
a cup suppressor bias in the range -300±45V. We esti- 
mated a ±0.8% beam flux integration uncertainty based 
on the difference of the good geometry (Faraday cup) and 
poor geometry (biased target arm) results. The beam 
was rapidly deflected from the target prior to and dur- 
ing arm movement. The timing cycle intervals were 
t x = t 3 = 1.50021 s, t 2 = 0.24003 s and i 4 =0.26004 s, 
and the (inverse) timing efficiency /3( 8 B) = 2.923 ± 0.005 
assuming ii/ 2 ( 8 B) = 770 ± 3 ms |J. 

In the limit of uniform beam flux, the 7 Be areal density 
is unimportant and the cross section is given by 



<j(E c 



Y a (E p )F a (E p )f3( s B) 
2N p N Bc (t)n/4ir 



(1) 



where E cm is discussed below, E p is the bombarding en- 
ergy, Y a (E p ) is the a yield above a threshold energy of 
895 keV, F a (E p ) is a correction for the fraction of the 
a-spectrum that lies below the threshold, N p is the inte- 
grated number of protons per cm 2 , NB e {t) is the number 
of 7 Be atoms and fi is the solid angle of the a-detector. 

In practice it is impossible to produce a completely 
uniform beam flux. To understand the error associated 
with this approximation, one needs to know both the 
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FIG. 1: Top panel: a-spectrum from 7 Be(p,7) 8 B at E cm = 
186 keV. Middle panel: 770 keV deuteron beam transmission 
ratios through different apertures, vs. raster amplitude. Bot- 
tom panel: 7 Li(d,p) 8 Li yield at 770 keV, normalized to the 
integrated beam flux through a 3 mm aperture, vs. raster am- 
plitude, measured with the same tune as the aperture ratio 
data. 

beam and target uniformities. It is particularly impor- 
tant that the target be confined within a small central 
area. This was insured by depositing the 7 Be on a Mo 
backing consisting of a 4 mm diameter raised post sur- 
rounded by a mask tightly pressed around the post, with 
post plus mask machined flat as one piece. After evapo- 
ration the mask was removed, eliminating unwanted tails 
on the 7 Be radial distribution [Q . 

The beam uniformity was determined by measuring 
the transmissions through 2, 3 and 4 mm apertures as 
functions of the (equal) amplitudes of the x and y trian- 
gular raster waveforms. Fig. [j] shows measurements with 
a 770 keV deuteron beam, and curves calculated by fold- 
ing a Gaussian with a rectangular function. The unifor- 
mity of the product of the beam and target densities was 
determined by the raster-amplitude dependence of the 
7 Li(d,p) 8 Li yield from the 7 Be target at E d = 770 keV, 
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FIG. 2: Top panel: 7 Be(a, 7 ) 11 C resonance profile corrected 
for small backgrounds due to cosmic rays and a nonresonant 
yield from 9 Be(a, n) 11 C. Bottom panel: 7 Be activity divided 
by the 7 Be decay curve, showing sputtering losses. 

shown in Fig. [l| The curve is a 1-parameter folding of the 
target density estimated from 7-activity scans, and beam 
profile determined by the transmission ratios, including a 
fitted target-aperture misalignment of 0.5 mm. The point 
at which this yield flattened out determined the mini- 
mum safe raster amplitude, and is similar to the point 
at which the aperture ratio data flattened out. We chose 
0.42 as the safe raster amplitude for 770 keV deuterons, 
and assigned a conservative ±1% nonuniformity uncer- 
tainty here. Aperture transmission curves, measured at 
most proton energies, determined the minimum raster 
amplitude for each energy and tunc for which the beam- 
target nonuniformity was < 1%. Independent estimates 
of the safe raster amplitudes were made by folding the 
target density distribution jl4| with beam-flux distribu- 
tions determined from the proton aperture-transmission 
data. 

-^Bc(i) was determined with the target arm vertical 
by counting 478 keV 7-rays in situ using a collimated 
Ge detector located on top of the target chamber. We 
assumed t 1/2 = 53.12 ± 0.07d |i| and a 10.52 ±0.06% 
branch |ll| to the 478 keV level. The Ge efficiency 647s 
was determined to ±1.3% from a fit to 14 lines from 
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125 Sb, 134 Cs, 133 Ba, 137 Cs and 54 Mn sources calibrated 
typically to ±0.8%(1ct) (l|, with jv = 2.2. We ob- 
tained a second 137 Cs source calibrated independently to 
The relative activity of the two 137 Cs 



±0.4%(1ct) 1 16 



sources agreed within ±0.1%. As can be seen in Fig. 
2.5 mCi of 7 Be was lost due to beam sputtering during 
the cross section measurements. 

TABLE I: Percent uncertainties AS17/S17. 



Statistical errors 


1.0-2.8 


Varying systematic errors: 




proton energy calibration 


0.2-0.6 


target thickness 


0.0-1.0 


target composition 


0.0-1.1 


Scale factor errors: 




beam-target inhomogeneity 


1.0 


integrated beam flux 


0.8 


target activity 


1.9 


solid angle 


1.2 


a-spectrum cutoff 


0.7 


backscattering 


0.5 


timing cycle 


0.2 


Total scale factor error 


2.7 



We inferred fl with the aid of a "far" Si detector 
located 47.42 ± 0.09 mm from the target and collimated 
to an area of 248.8 ± 0.4 mm 2 . From geometry, £lf ar — 
0.1078 ± 0.0004 sr, where the zero of the distance scale 
was checked using a 148 Gd a-source. Sl/Qf ar was de- 
termined using the 7 Li(d,p) 8 Li reaction. A differential 
correction for a-particles lost below the threshold was 
applied based on the (d,p) angular distribution [jlTj and 
SRIM Q calculations including 8 Li-straggling. We ob- 
tained = 3.82 ± 0.04 sr. This result was checked using 
different detectors and different size collimators for f2/ r- 

The yields Y a (E p ) were corrected for a small beam- 
off background (3.9% at the lowest E p ). The beam- 
related background was checked at several energies and 
found to be negligible. The cv-spectrum cutoff factors for 
7 Be(p,7) 8 B were estimated from SRIM calculations, in- 
cluding straggling, fitted to 23 different spectra. F a (E p ) 
varied linearly from 1.039 ± 0.007 at E p = 221 keV to 
1.086 ± 0.008 at 1379 keV. The accelerator energy cali- 
bration was determined to ±0.17% from 19 F(p,ev7) 16 
resonances at E p = 340.46 ± 0.04, 483.91 ±0.10 and 
872.11 ± 0.20 keV @. 

Corrections for energy averaging of the proton beam 
due to finite target thickness are important, particularly 
at low E p . We directly measured the beam energy loss 
profile in the target using the narrow (r << 1 keV) 
7 Be(a, 7) n C resonance [^0| which we found at E Q = 1378 
± 3 keV. The mean a-energy loss was 26 ± 2 keV, based 
on the average of three measurements, one of which is 
shown in Fig. |^. The excellent reproducibility of the ap- 



parent 7 Be(a,7) 11 C resonance energy measured in the 
middle of, and after the 7 Be(p,7) 8 B measurements (AE a 
= 1 ± 3 keV), indicated negligible carbon buildup and 
target damage due to bombardment. 

An important error in some previous experiments was 
loss of 8 B from the target due to backscattering (and loss 
of 8 Li when 7 Li(d,p) 8 Li was used for absolute cross sec- 
tion normalization) [^l], |22| . These losses may be sizeable 
when a high-Z backing is used, or if there are high-Z con- 
taminants in the target. We made the first direct mea- 
surements of 8 B backscattering losses in the 7 Be(p,7) 8 B 
reaction using our 7 Be target in a fixed mount, and large- 
diameter water-cooled Cu catcher plates on each end of 
the rotating arm. A 4 mm hole in the center of each 
plate allowed the beam to pass through. We found small 
backscattering losses of 1.3 ± 0.3% and 0.9 ± 0.2% at E p 
= 724 and 1379 keV, respectively, and made a constant 
1.0 ± 0.5% correction to our data for this effect. 
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FIG. 3: Si7(-B cm ) vs. Ecm from this work. Error bars are 
statistical plus varying systematic errors. Solid curve: DB 
theory plus a Breit-Wigner resonance. Dashed curve: DB 
theory. Inset: resonance region. 

Fig. U shows our S-factors calculated from the relation 
S 17 (E cm ) = a(^ cm ) J B cm exp[(S G /S cm ) 1 /2] ( see e .g. §) 

with Eg = 13799.3 keV. We computed E cm by invert- 
ing the expression a{E cm ) — d, where a was obtained 
by fitting the cross section data including averaging over 
the target profile, and a is the corresponding unaveraged 
cross section. These E cm values are very close to the 
mean proton energy in the target, except near the res- 
onance where they differ by < 1%. Fig. [3] also shows a 
fit to all our data of the (scaled) cluster model theory 
of Descouvemont and Baye J23j plus an E cm = 630 ± 2 
keV Breit-Wigner resonance (with energy-dependent T p 
and T 7 ). This fit yields Si 7 (Q) = 22.5 ± 0.6 eV-b and 
X 2 /^ = 1-3 (^=25) p3, where the quoted uncertainty 
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includes the scale factor error of ± 2.7 % (Table |). Fits 
with other theories p5| did not reproduce our measured 
energy dependence as well (x 2 /^ = 1.7-16). 

The theoretical uncertainty in the energy dependence 
of S17 decreases with beam energy below the resonance, 
as the capture becomes increasingly extranuclear. There- 
fore it is important to determine Si7(0) from low energy 
data. Fitting the DB theory to our data at E cm < 300 
keV we find Si 7 (0) = 22.3 ± 0.7 eV-b and x > 2 )v = 0.3 
(here, as above, the error includes statistical plus system- 
atic contributions) . In addition, there is an extrapolation 
uncertainty, which has been estimated to be as small as 
±0.2 eV-b ^BJ, and which we estimate conservatively as 
± 0.5 eV-b from the rms deviation of 11 different theo- 
retical fits to our data for E cm < 300 keV @ . Thus our 
final result is 

Si 7 (0) = 22.3±0.7(expt)±0.5(theor) eV-b. (2) 



30 - 
28 



CO 



21 



Parker □ Hammache A 
Kavanagh o Strieder ♦ 
Filippone A this work ▼ 





i i 

-L i 


— ! : 

• 





FIG. 4: Si 7 (0) from our fits of the DB theory to E cm < 425 
keV data from this and previous measurements. Horizontal 
lines indicate the 19 t| eV-b range recommended by ^8] . Fits 
over a wider i? p -range give similar results but with smaller 
errors for other experiments. 

In order to compare all direct measurements below the 
resonance, we made DB fits to all data at E cm < 425 
keV - this work and ||, ||, 0, p, [H], renormalized to 
cr[ 7 Li(d,p) 8 Li] = 152 ± 6 mb p9] where appropriate. The 
results are shown in Fig. ^. Results from |3|, [|, |?J may 
suffer additional error from 8 B and 8 Li backscattering 
losses; in B, calculated corrections were applied, while 
in 0, a low-Z backing was used and losses were assumed 
negligible. 

We have reduced the error on Si 7 (0) so that it no longer 
dominates the uncertainty in the calculated solar 8 B v e 
production rate. While our Si 7 (0) value agrees within 
errors with the previously recommended value of 19 + 2 
eV-b ^8|, it is 17% larger. Thus 17% more of the 8 B solar 
i/ e 's oscillate into other species than given in ref. B. 
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